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Abstract 

We extend the method discovered by A Y Alekseev et al [|l| to the case of 
fermions in external fields. A general formula for conductance G is proved. In 
the (1+1)-D case with symmetry at time reflection, it is shown that: G = e^ /h 
+o (a^j, where a is the strength of the external field. In (3+l)-D free case, 
it is checked that G = ne^ /h, where n is the number of the filled energetic 
bands of the transversal quantization. 

I. INTRODUCTION 

It was proven relative recently that the equilibrium statistical mechanics may be used at the 
description of some transport phenomena in solids. Using this approach, A Y Alekseev et al 
]![] have found a compact formula for conductance. We apply in this paper their method to 
study the transport properties of the fermionic systems. The method have been applied to 
the cases where the backscattering processes were not present. In this paper we show how 
to incorporate in the theory the backscattering processes. The presence of impurities can 
be modeled by introducing an external field. To analyze fermionic systems in external fields 
we make use of the implementation of the one particle picture. This will enforce us to work 
only in the fermionic representation. 

We calculate first the conductance of the free (1+1)-D fermionic system and rediscover 
the results of [Q, obtained by using the bosonic representation. Then we find a natural 
way to extend the formalism to the case of the fermionic systems in external fields. A 
general formula for conductance will be deduced. It is shown that the first correction to 
the conductance is quadratic in the strength of the external field. This is consistent with 
the fact that the conductance should be lower than e^ /h for both cases: attractive and 
repelling external fields, because the backscattering processes are present in both cases. 
However, because of the complicated analytical expressions, we are not able to check that 
indeed the conductance is strictly lower than e^ /h. We anticipate interesting behavior of 
the conductance with the increase of the external field strength: discontinuities when bound 
states appear or disappear. 

The last section treats the (3+l)-D fermionic systems. It was argued in the literature 
0] that the conductance of a pure quantum wire is proportional with the number of filled 
energy bands of transversal quantization. We check this for the free (3+l)-D fermionic 
system. For (3+l)-D fermionic systems with external interactions only in the directions 
perpendicular to the measured current, one can combine our results with those of [^ to 
study the quantization of the conductance for Quantum Hall effect. 



II. THE GENERAL PRESCRIPTION 

We will briefly discuss first the method discovered by A Y Alekseev et al in [|I|. Given a 
quantum theory described by a Hamiltonian H, the first step is the identification of two 
physical components of the quantum field, the left and right parts, and the corresponding 
charges: Ql,r- They must satisfy the following commutation relations: 

[H,gL,K]=Oand [Ql.Qr] = Q. (1) 

Denoting the Noether conserved currents corresponding to the charges Ql,b. by Jl,r and the 
total current by j, the expectation value of the current, in the equilibrium state correspond- 
ing to the chemical potentials hl^r (denoted by (•)„), is: 

(J iy)), = -^ {^^L - f^R) I dx { [j° (x) , a (y)] >^ . (2) 



Here, a is the solution of: 



f = eVa , j = edta, (3) 



and ja denotes the axial current, Jl — Jr- The electric potential can be expressed as: 

eV = ^L- fJ-R- 

In the following, we will refer to the above method as the general prescription. If one 
works in the fermionic representation, the commutator [pa [x) , a {y)] cannot be calculated 
directly and we apply the following strategy. We calculate first the expectation value: 
f {x,y) = {[pa{x) , p(l/)])„- Then we search for a vector F{x,y) such that: \/yF{x,y) = 

f (x, y). In the case of lim F [x, y) exists (and is well defined in the sense that it does not 

depend on the direction), the expectation value of the commutator [pa{x),a{y)] is given 
by F{x,y) — lim F{x,y). This assures that pa{x) and a(?/) decuple when x and y are 

separated by infinite large distance, which must be the case if we assume that a is a local 
observable. Along this paper, we will keep the notation / and F for the above quantities. 



III. (1+1)-D FERMIONIC SYSTEMS 
A. (1-|-1)-D free fermionic system 

Let us consider the free fermionic system described by the Hamiltonian: 

k 

for some finite length, L, and periodic boundary conditions. Following the general prescrip- 
tion, the first step is to define the left and right components of the field. We propose the 
zero time left and right fields and the corresponding charge operators to be: 



--^XnuJC 



^L,R {X) = J2 -Tf'^k , Ql,R = J2'^l ^'^- (5) 

It is easily to check that the following commutation relations held: 

[Ho,gL,ij] = o, [QL,g«] = o. (6) 

So, we are in the conditions required by the general prescription. The axial charge density, 
Pa-, is defined by: 



Ql-Qr= dx {pL (x) - PR (x)) = dxpa (x) . (7) 

The left and right charge densities are not uniquely defined. We choose the following repre- 
sentation: 

PL,R (x) = \ [^^ {TIl,rx) ^ {x) + ¥ {x) ^ {IIl,rx)\ (8) 

where YiL,R represent the projectors on the left states: /c > 0, respectively on the right states: 
A; < 0. The notation is coming from: 

pikx n pik(x—y) 

^L,«(:r) = E ^f'^k = E / dy——^{y) 

= JdyUL^ji{y,x)^iy) = ^iq), (9) 

where q is given by: q (y) = [IIl,rS {x — ■)] (y) and we have denoted: IIl,rS {x — ■) = IIl,rx. 
The zero temperature expectation values corresponding to the chemical potentials pi r = P^ 

eV 

— , can be calculated by employing the functional integral over the Grassmanian variables 

with the measure corresponding to the covariance: 

Ckk' = hk'-, rryTT^ ttt tttT) (10) 

iko - [k'^/2m - PL X> [k) - Pr X< {k)\ 

where x^ ^-^e the characteristic functions of (— oo, 0] and [0, oo) respectively. It follows from 
Appendix A that the expectation value of the following commutator is given by: 

{[pa (x) , p (y)])^ = tlm[{^ {{Ul - Ur) x) , ^>^ {y) ] {^>^ {x) ^ (y)) J (11) 

and further: 

(k {x) , p (y)])^ = I im [(n^ - n^) (y, x) {yIl ^[o,^.] + ^r ^[o,^«]) {x. y)] ■ (12) 

We have denoted by E^ the spectral projection of Ap (Ap the Laplace operator on L 

with periodic boundary conditions) for some Q E R. In the infinite limit {L -^ oo) , we have: 



-i sin ki (x — y) + sin k^ (x — y) 
2v^ (x - yf 



f (x, y) = ^^ -^-'^L,^-yj^^^^^^r.n,^-yj _ ^^3^ 



Following our strategy, we find: 

ZTT^ l^ X — y 

+ ^ j/c.? Ci [A:^ (x - y)] fh^j ' (1^) 

with the trivial limits lini F {x, y) = 0. From formula H, the current density is: 

if 6 

(i (2/))/. = -^ if^L - f^R) dxF {x, y) = [fiL - fin) - (15) 

which leads to G = I/V = e'^/h. The main goal of this subsection was to check the validity 
of our strategy. 



B. (1+1)-D fermionic systems with external interactions 

We will consider in this subsection (1+1)-D fermionic systems in external fields. We restrict 
here at the second quantization of the one particle problems described by the Hamiltonians of 

the following type: H = A+V, considered again on a finite length, L, and with periodic 

2m 
boundary conditions. We will consider first a compact supported potential, suppV C [a,b], 

with —L/2 < a < b < L/2. Also, we consider the case when the Hamiltonian H doesn't 

have any bounded state in the limit L -^ 00. In the above conditions, the eigenvalues of the 

Hamiltonian are again — ^, with kn in a discrete set of numbers. Let e^ be the eigenvector 

2m 

corresponding to some k with e^ (x) oc e"*'^^ for x < a. A; > and e^ (x) oc e*'^^ for x > b, 
k < (the "scattered" states from left and right 0). The second quantized Hamiltonian is: 

k 

where "^k has now the following meaning: 

^fc= / t/xefc(x)^(x). (17) 

The left and right components of the field and the corresponding charge operators can be 
defined analogous: 

^L,/? (x) = J2el (x) vP, = v]/ (Ul,rx) , Ql,r = J2 ^l^k, (18) 



where IIl,r are again the projectors on the states with k > and A; < respectively. It is 
straightforward that the charge operators commute with the Hamiltonian. So we are again 
in the conditions of the general prescription. The current density will be given by the same 
expression |^. Further, because {"^k} satisfy the same anti-commutation relations we can 
conclude that: 



{[pa (x) , p (y)]) , 



(19) 



2 Im [{vl/ ((n^ - n^) x) , vl/t (y)] (vpt (a;) ^ (^))^ 
According to Appendix B, the above expression can be written completely analogous as: 

{[pa (x) , p (y)]) = I Im [{Hl - TIr) {y, x) {Ul ^[o,;.,] + Hr Ep,^^,) (x, y)] , (20) 

where Eq is the spectral projector of H corresponding to some Q & R. In the infinite limit 
(L -^ oo), we can use the results of the last section if a supplementary condition is fulfilled, 
which is: the Moller operators W± = s — lim q^th^-iTHq g^ist and are complete. In this 

case one can use the interwining properties of Moller operators, Eq = W+ Eq Wl. and the 
obvious relation: II^^jj = PV+n)^^iy|. In more complicate situations, the last relation can 
be used as the definition of the left and right states. The index (0) is referring to the free 
case. Using the explicit expressions of the projectors for the free case, one has: 



{[pa{x),p{y)])^ 



ilm 



W^AW. 



{y,x) 



W^BW^ 



{x,y) 



where A and B are the operators corresponding to the kernels: 



A{x,y) 



1 



iTi {x — y) ' 



B{x,y) 



sm 



kL+kji 



X 



y)] 



'■-{x-y) 



TT{x-y) 



(21) 



(22) 



One can go one step further by employing the following argument. If one changes the sign of 
the applied electric voltage, then the new current must have the same intensity but opposite 
direction. This is due to symmetry at time reflection. It implies that (j (?/))„ must be anti- 
symmetric at the interchange: k^ ^ k^. Because in the formula |^ there is already the term 
{Pl — Pr) in front of the integral, it follows that we should retain in {[pa {x) ,p (y)])„ only 
the symmetric term relativ to the interchange: ki ^ kj^. From the expression of B {x,y) it 
follows: 



{[pa (x) , p (y)]) , 



ilm 



W^AWl 



{y,x) 



W^ ReBWl 



ix,y) 



(23) 



We consider that the formula (|23|) can be used to calculate the conductance for more general 
cases (for example when bounded states are present). Indeed, the transport of electric 
charges can be done only by scattering states. One can see that the formula (^) takes into 
account only the scattering states so it can be used directly without going back and redefine 
the currents (which should incorporate only the scattering states). 

C. The dependence on the strength of the external field 

Let us introduce a scaling factor, a, in front of the potential: V —>■ aV. We study in this 
subsection the dependence of the conductance on the parameter a. The principal result is: 
dG/da\^^Q = which leads to the following expansion of G: 



6 



G=- + G^'^a' + o (a') 
h 



(24) 



One could anticipate this by considering both cases, ±V. If the first perturbative correction 
were non-zero, then, for one of these two cases, the conductance were greater than that of 
the free case which is wrong because in both cases the backscattering processes are present. 
This argument also implies that G^^' should be a negative number. 
The stationary representation of the MoUer operator is P]: 



W^ 



a 



lim / i?o (A - is) F i?o (A + is) d\ 

2711 £\0 



a 



2m 



oo 

lim / /?o (A - ie) VRa{X- ie) V {Rq (A + ie) - i?o (A - ie)} dX 

e\Oj 



= 1 + aWi + a^W2, (25) 

where R^ {z) represents the resolvent of —A + aV . Then one can see that: 
d 



da 



[pa {x) , p {y)] 



zim 



a=0 



{[ 



Wi,A 



{y,x) B{x,y) + A{y,x) 



Wi,B 



{x,y)]. (26) 



With the observation of the last subsection: 



^[Pa{x),p{y)] 



a=0 



ReWi, A 



{y, x) Re [B {x, y)] + A (y, x) 



ReWi, ReB 



{x,y) 
(27) 



In our further calculations, we will deduce for Wi a formula similar to that of [Q. The 
calculations are lengthy and are presented in the Appendix C, where it is proved that 

-T- [Pa {x) , p (y)] = 0. 
da „^o 

The coefficients of the upper powers of a are to complicated (at least for the author) to 

be handled analytically. However, a numerical approach will not be a trivial thing and it 

can make alone the subject of a new study. 



D. Comments 



We have pointed out before, that the formula ^ can be used in more general context. The 
Moller operators project out the bounded statets (and zero energy resonances). As a result, 
only the scattering states contribute to the current. One can see that there is a major 
difference between the cases with and without bounded states. We expect discontinuities 
of the conductance (versus strength of the external field) when bounded states are forming. 
Thus, the investigation of the transport properties may be a complementary way to the 
method of scatterings in detection and measuring of the spectrum of energy for quantum 
systems. 



IV. (3+l)-D FERMIONIC SYSTEMS 

The primary goal of this section is to prove that the conductance of the free (3+l)-D 
fermionic system is proportional with the number of filled energy bands of the transversal 
quantization. The left and right fields and corresponding charge operators can be chosen as: 

^L,H (^) = ^ E ^^'"^A? , Ql,r = E ^ ^.^, (28) 

where Xi denotes the direction of the measured current and v is the volume. One can follow 



the same steps as in (1+1)-D case to prove the validity of the formula |TT]. According to 
Appendix D, in the infinite limite: 

f{x,y) = '^^''^~y'-\ kiu{kUxi-yi)) + k''^u{kn{x,-y,))}, (29) 

71 

, , , cosx sinx ^ , ^ , . ,. . , 

where u [xj = — — . Because the current nows only m xi direction, we can choose 

the vector a along xi and we can solve easily the equation \7yF = f (x, y): 

F,{x,y) = {[pa{x),a,{y)])^= '^^''^~y^^ {klUikL{x^-y^)) + klU{kn{x^-y,))}, 

(30) 

• 1 T T / \ COS X vjl X Sm X Sm X tt. i roOrr/\7 /r\* rn 

with U [xj = 1 5- H . Using that /_ U fx) ax = vr/z, it follows: 

3x^ 3 3x'^ 3x 

(Ji (2/))^ = {[Ql- Qr, ai (y)])^ = / ^"^^ / ^"^1 ^^P- (^) ' ^ ^^)]^ = ^ (^' + ^^) = '^^ ^^^) 

which leads to: 

f if f 

I = J dy^{3^{y)) = ~—{iiL-liR)ikl^S = -{nkl/\S)-V. (32) 

Thus: G = I/V = (TrkpAS) —. Of course, the pre-factor nkpAS represents the number 

of the transversal conducting channels. The result may be generalized to the case when 
transversal interactions are included. In [Q it was calculated the change in the number of 
energetic bands below the Fermi energy for 2-D quantum systems, when infinitesimally thin 
magnetic flux tubes are introduced. We think that the two results can be combined to study 
the quantization of the conductance for Quantum Hall effect. 
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V. APPENDIX A 



In this section we will calculate explicitly the commutator for the free fermions: 
[pa {x) , p {y)] = i [M/t ((n^ - n^) x) vj/ (x) , vl/t (^) xi/ 



+ - [¥ (x) ^ {{Ul - Ur) x) , ¥ (y) ^ (y)] . (33) 

We can use here the general formula: 

[ab , cd] = a[{b , c} d — c{b , d}] + [{a , c} d — c{a , d}] b (34) 

which leads to: 

[*t ((n^ _ n^) a;) ^ (x) , vj/t (y) ^ (y)] = S {x - y) ^^ ((H^ - U^) x) ^ (y) 



-{vl/t((ni-n^)x),vl/(y)}il/t(^)xl/ 



X) 



and 



(35) 



[vl/t (x) ^ {{Ul - Ur) x) , v&t {y) H/ {y)] = {M/ ((n^ - n^) x) , M/t (y)} M/t (x) ^ (y) 



-6{x-y)^Hy)'^{{^L-TlR)x). 
In consequence, the commutator can be expressed as: 

[Pa (x) , p {y)] = \6{x- y) [M/t ((n^ - n^) x) ^> {y) - M/t (y) M/ ((n^. - Hr) x)] 



(36) 



+ ^ {^ {{Ul - Ur) x) , Vj/t {y)} H/t (x) H/ {y) 



- - {M/t ((n^ - Ur) x) ,^{y)} ¥ (y) ^ (x) . 
By a simple operation, one can put the first term in the form: 

5{x-y) [^[ (x) ^i? (x) - ^^ (x) ^L (x) 



(37) 



(38) 



and one can see that its expectation value is zero. If now one considers the expectation 
value of the next two terms, the result can be written in the compact form: 



{[Pa (x) , p (y)])^ = z Im {vl/ ((n^ - Ur) x) , M/t [y)] (M/t (x) M/ {y)) ^ 



(39) 



by taking into account that the anti-commutator is just a number. Moreover, by using the 
field anti-commutation relations: 

{^ ((n^ - n^) x) , vj/t (j/)} = 1 ^ e^^-y^ _ 1 ^ e^M-.) = (n^ _ n^) (^, ^) (40) 

fc>0 fc<0 

which in the continuous limit becomes: 

1 roo -1 ^00 ^00 ikrj i 

-jjk (,> (.) - x< (*))e--" ^ ^J_Jt lj„-e'^-'^ ^ -^-. ,41) 
The propagator can be calculated by using the functional integral: 

(vl/t (x) vl/ (y))^ = y vl/t (a;) vl/ (^) rf^^ (42) 

with Grassmann Gaussian measure corresponding to covariance: 

^'^ iko - (A;V2m - fn X> W - [Xr x< {k)) ' 



(vl/t (x) M/ (y))^ = ^ [k^/2m - /Xi x> (k) - m X< (k)) 



^-ik(x-y) 



= (nL^[o,^^] + ^rE[o,i,^]) [x, y) . (44) 

In the infinite limit, the propagator becomes: 

kL 

(vl/t (3;) vl/ (^)) = /■ ^ e-./c(.-.) ^ sinfcHa:-!/) .Ji^^(.-,)_ ^^ 

^ ^ ' ^^"^' J 27T 7r(x-?/) ^ ^ 

-kn 



VI. APPENDIX B 

We prove here that the formula |T2| is still valid for the interacting case. Because the fields 
^k satisfy the usual commutation relations, < ^k' ; ^I [ = Sk,k', the formula ^ results by 
following the same steps. Now 

{^iiUL-nn)x),¥iy)}= J2 ekiy)el,ix){^k',^l}- Yl ek (y) e*,, (x) {^ k' , ^l] 

k,k'>0 k,k'<0 

= E ^'^ (^) ^fc (^) - E ^'^ (^) ^fc (^) = (^L - n^) [y, x) . (46) 

fc>0 fc<0 

For the expectation value (\l/^ (x) \l/ (|/)) can be calculated in the same manner as in Ap- 
pendix A. The covariance of the Grassmann Gaussian measure for the interacting case has 

10 



the same matrix elements but relative to the orthonormal basis {cfc}^. The functional inte- 
gral leads to: 

{¥ (x) ^ (y))^ = 5^ e {e/2m - fXL X> (k) - fxn X< (k)) e^ (x) el (y) 



(Hl^ICa^l] + ^REio,m]) i^^y) ■ 



(47) 



We can conclude that the expression 12 is still valid for this case. 



VII. APPENDIX C 



dG 

We show here that — — 

da 



0. 



a=0 



Wi (x, y) 



f27r) 



OO /"OO POO 

d\ dz dk 



— OO ^ — OO ^— OO 



-ik{x—z) 



A;2 - (A - ie[ 



-Viz) 



OO ik'(z-y) 

dk'- 



k'^ - (A + ie) 



CO /"OO 



(27r)"^/^ / dA; / rffc 



— OO >/ — OO 



, g-^'^^l/ (fc' - k) e'^^'y 
fc'2 - fc2 - 2ie ' 



(48) 



where indicates the Fourier transform. We make the following change of variables: 



i = k'-k 
f] = k' + k 



(49) 



d (k k') 
with the Jacobian: -rrr^ — r = 1/2. Then: 
5(e,^) 



W.{x,y) = -^^^l dul d,,'- ^^^ 



OO J — OO 



:27r) 



-3/2 



°o Qi£.{x+y)/2Y /^\ /-OO ^_j^(a;_y)/2 



dC 



e 



(ir/ 



r/ - 2ie/^ 



i{27T 



,-1/2 



c?^ [X> {x - y) x> ii) -x<{x- y) X< (0] 



,iax+y)/2Y (^^^ 



(50) 



which is the expression of Wi in [||. Further: 



ReVri(x,y) 



i(2n] 



-1/2 



c?^ [X> (a; - y) X> (0 - X< (a; - v) X< (0] x 
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X- 



2^ 



i(27i) 



-1/2 



[X> {x-y)- X< {x - y)] / d^ x> (0 



^i^{x+y)/2Y ^^^ _ ^-i^{x+y)/2y (^_^^ 

2e 



(2r^-^/^ 



ill 



oo ^i({x+y)/2Y /^\ 

[X> (x-y)- X< (x - y)] I d^ - 



We write generically: 

Re Wi {x, y) = [x> {x - y) - x< (x - y)] Q {x + y) 



(51) 



(52) 



The general expression for the commutators ReWi,C , where C{x,y) = C{x — y), will 
be: 



Rel^i,C 



(x, y) = dz [x> {x- z)- x< {x - z)] Q{x + z)C{z-y) 



dzC{x- z) [x> {z-y)- x< {z - y)] Q {z + y) 



dz [x> {x - y - z) - x<{x - y - z)]{Q {x + y + z) - Q {x + y - z)) C {z) (53) 



after a change of variable: z ^ z — y in the first integral and: z — > —z + x in the second 
integral. Concretely: 



ReW^i,C 
It follows: 



[x,y) 



{2n) 



-1/2 f°° , e" 



({x+y)/2y (c\ / rx-y /-oo \ /^^ 

dC- ^^^(/^ -[ )dzsin[^%^]C{z). (54) 



oo J x—y 



KeWi.A 



ix,y) 



(2^) 



-1/2 f^ . e'' 



j^^^"-'>'^-vio (r_r), ;"'y^2 



f 



— oo J x—y 



l-RZ 



i(2tt) 



Tj2 



oo ^ii,{x+y)/2y /^\ 

d^ [Si {x — y) — Si (— oo) — Si (oo) + Si (x — y)] 



.,2.wi J^ 1 si(-(x-y) 



(55) 



and 
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ReWi,ReB 



ix,y) 



(2n)-y- r ^^ e'^^^+yy^vjO 



"-^ /"" ^ . f^z\ sin (kLz) + sin (km) 
az sm I — I 

—CO J x—y ' \ ^ I ZTIZ 



2(27r) 



^^ e-€(^+^)/V (^) ^ |-s/ _ f^^^^ cos [(fci - e/2) ^] - cos [{kr. + i/2) z] 



oo J x—y 



2z 



+^^ /" rf^ e'^^-+^^/ V (0 1^ /"-^ - Ddz '°' f^^^ ~ ^/^^ ^] " '°' f^^^ + ^/^^ '^ 



2(27r)^/2 



e 



oo J x~y 



2z 



2(27r) 



^ / ^^ ^^^ {Ci [(^L - e/2) (x - y)] - Ci [(fc^ + e/2) (x - y)] 



oo 



i 



+ Ci [(fc« - i/2) [x - y)] - Ci [{kn + e/2) (x - y)]} . 
Now we start to apply the strategy of section 2. We need to calculate: 



(56) 



dx / dy 



W^,A 



{y, x) Re [B (x, y)] + A (y, x) Re M/i , 5 (x, y) . (57) 



Changing the order of integration, one can see that the situation is similar to 



dxe^^"" I dye 



<^^-v)l^g (x - y) 



dx V i^ 



dxe-'^^^-y'>/^g{x-y) 






dxe-'^^^-y^/^g{x-y) 



a;=-oo ^— oo 



oo gj^x 



+ / dx — e-*«(^-^)/2^ (x - y) . (58) 



The first term is zero for our concrete case. The integral may be calculated explicitly [Q. 
Because the current is conserved, we can take any particular value of y. Choosing y = the 

above integral becomes: /^ dx — . Applying this scheme to 0, the first correction 



to the current density is: 



^^ 



le 



^•^'^'^>A^ = ~i^^^~^^M^(^ 



^, no /"^ ^^ e-g"/^[sin(A;^x) + sin(A:^x)] g. ( ix 



^e 



2ttx 



--^ r d^^ r dx{Ci[{kL-^/2)x]-Ci[{kL + ^/2)x]} 

^ ' J —oo ^S J— oo 



^iix/2 



mx 
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+^dw^ ^^^ dx{Gi[{kn-i/2)x]-Gi[{kn + i/2)x]}^). (59) 

We concentrate now at: 

X \2 ./__ X V 2 



^1 ^^ cos ((fcz. - e/2) x) - cos ((fcz. + e/2) a;) g. /^a: 

, AX \ A 



^°° rfx ^ {Ci p^ - e/2) x] - Ci p^ + e/2) x]} Si (^ 



2 



^{Cip,-e/2)a:]-Ci[(A:, + e/2)x]}Si[^ 



sm 



- - / t/x {Ci [{kL - e/2) x] - Ci [{kL + e/2) x]} ^^. (60) 

Finally: 

/ rfx^ siin^z^g.^^ 1 / ^^e |ciPz.-e/2)x]-CiPi + e/2)x]} 

J —oo •^ \ '-' / ^ J —oo •^ 

(61) 

and one can see now that in ^ terms are canceling one each other, the final result being 
zero. During our calculations, we have changed the order of integration by applying the 
Fubini's theorem. As is pointed out in 0, we have to be carefully especially with oscillatory 
integrands. With a little effort, one can show that we are in the case where the Fubini's 
theorem can be applied. 

VIII. APPENDIX D 

In this section we will calculate the formulas used in section IV. By using the field 
anti-commutation relations it follows: 



{ vp ((n^ - n^) x) , vl/t (^) } = 1 ^ e^H^-y^ - i ^ e^^(^-^) = (Hz. - Ur) {y, x) (62) 
and in the infinite limit: 



V ^^ V 

fci>0 fci<0 
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{^ ((n^ - n^)a:) , vl/t (^)} = ^^jd^k (x> (A;i) - x< (h)) e^'^^-y^ 



(27r) 

= —f \^ {x^ - va) , (63) 

m [xi - yi) 

where we have denoted by x\_ the coordinates: X2 and 3:3. Also, 

1 Lc" i ^r* -ft \ 

{¥ {x) M/ (y))^ = ^ {e/2m - fiL X> ih) - f^R X< ih)) ^ 



k 

= {liLE^,^,,^ + IiRE^,,,^{) {x,y) . (64) 

In the infinite hmit, 

— e {e/2m - /i^x> ih) - f^RX< ih)) e-^^(^-^). (65) 

In consequence: 

{[Pa (x) , p (y)]), = ^ Im ~^ ["^^ ' ^^^ / ^ e (A;V2m - /.^ x> ih) - ^,n X< ih)) e~^'^^-y^ 
^ tiT {xi -yi) J 2-K ^ 

iS {x± — y±) f (Pk 



TT {xi - yi) J 27r 



/(rk 
— [k'^/2m - fiL X> (h) - m X< (h)) cos ki {xi - yi] 



i5 {x± - y_i_) [^ dk 



fcL rii. _ /-Vs 



7r(^xi -yi 
iS {xj_ - yi_) /■''^ dk 



/ — k'^2TT dO sin 9 cos {k(xi — yi) cost 

Jo 27r Jo 



f " dk r 

I — k'^2Tc d6 sin 6 cos {k {xi — yi) cos t 

Jo 27r J^/2 



TT {xi - yi) 

_ iS ix± - y±) \ ki cos ki {xi - yi) + kn cos kn {xi - yi] 
Ti{xi-yif I xi-yi 

sin kL {xi -yi) + sin kn {xi - yi] 



{xi - yif 



We conclude: 



f {x.y) = —{klu{kL{xi -yi)) + k\u{kR{xi - yi))] , (66) 

TT 



cosx smx 
where u [x) 



x'^ x^ 
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